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ABSTRACT 

We use absolutely calibrated data from the ARCADE 2 flight in July 2006 to model Galactic 
emission at frequencies 3, 8, and 10 GHz. The spatial structure in the data is consistent with a 
superposition of free-free and synchrotron emission. Emission with spatial morphology traced by the 
Haslam 408 MHz survey has spectral index /Jsynch = —2.5 ±0.1, with free-free emission contributing 
0.10 ± 0.01 of the total Galactic plane emission in the lowest ARCADE 2 band at 3.15 GHz. We 
estimate the total Galactic emission toward the polar caps using either a simple plane-parallel model 
with esc |6| dependence or a model of high-latitude radio emission traced by the COBE/FIRAS map of 
Cil emission. Both methods are consistent with a single power-law over the frequency range 22 MHz to 
10 GHz, with total Galactic emission towards the north polar cap Tba! = 10.12 ± 0.90 K and spectral 
index 0 = —2.55 ± 0.03 at reference frequency 0.31 GHz. Emission associated with the plane-parallel 
structure accounts for only 30% of the observed high-latitude sky temperature, with the residual in 
either a Galactic halo or an isotropic extragalactic background. The well calibrated ARCADE 2 maps 
provide a new test for spinning dust emission, based on the integrated intensity of emission from the 
Galactic plane instead of cross-correlations with the thermal dust spatial morphology. The Galactic 
plane intensity measured by ARCADE 2 is fainter than predicted by models without spinning dust, 
and is consistent with spinning dust contributing 0.4 ±0.1 of the Galactic plane emission at 23 GHz. 
Subject headings: radio continuum: ISM, radiation mechanisms: non-thermal, cosmic microwave 
background 


1. INTRODUCTION 

The cosmic microwave background (CMB) is a valu- 
able probe of physical conditions in the early universe. 
Its frequency spectrum records the history of energy 
transfer between the evolving matter and radiation fields 
to constrain the energetics of the early universe. We view 
the CMB through diffuse emission from the interstel- 
lar medium. At centimeter wavelengths, the dominant 
contributions are from synchrotron emission originating 
from electrons accelerated in the Galactic magnetic field, 
and free-free emission (thermal bremsstrahlung) from 
electron-ion collisions.. Thermal emission from interstel- 
lar dust is negligible, but electric dipole emission from a 
population of small, rapidly rotating dust grains could 
contribute a substantial fraction of the total Galactic 
emission at wavelengths near 1 cm (Draine & Lazarian 
1998; Dobler & Finkbeiner 2008; Mivllle-Deschenes et 
al. 2008). Absolute measurements of the sky tempera- 
ture at centimeter wavelengths measure diffuse emission 
to separate Galactic from primordial emission and pro- 
vide information on physical processes in the interstellar 
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medium. 

The Absolute Radiometer for Cosmology, Astro- 
physics, and Diffuse Emission (ARCADE) is an instru- 
ment to measure the absolute temperature of the sky in 
search of distortions from a blackbody spectrum. AR- 
CADE operates at centimeter wavelengths between full- 
sky surveys at radio frequencies below 3 GHz arid the Far 
Infrared Absolute Spectrophotometer (FIRAS) survey at 
frequencies above 60 GHz. It consists of a set of cryogenic 
radiometers observing at 37 km altitude from a balloon 
payload. Each radiometer uses a double-nulled design, 
measuring the temperature difference between a corru- 
gated horn antenna and an internal reference as the an- 
tenna alternately views the sky or a full-aperture black- 
body calibrator. The internal reference can be adjusted 
to null the antenna-reference signal difference, while the 
calibrator temperature can be independently adjusted 
to null the sky-calibrator signal difference. ARCADE 
thus measures small shifts about a precise blackbody, 
greatly reducing dependence on instrument calibration 
and stability. The calibrator, antennas, internal refer- 
ence, and radiometer front-end amplifiers are mounted 
within a large liquid helium Dewar and are maintained 
near thermal equilibrium with the CMB. Boiloff helium 
vapor, vented through the aperture, forms a barrier be- 
tween the instrument and the atmosphere to allow op- 
eration in full cryogenic mode with no windows between 
the optics and the sky. 

A 2-channel prototype (ARCADE 1) flew in 2001 and 
again in 2003, observing the sky at 10 and 30 GHz to 
demonstrate the feasibility of cryogenic open-aperture 
optics (Kogut et al. 2004; Fixsen et al. 2004). A second- 
generation 6-channel instrument (ARCADE 2) flew in 
2005 and again in 2006. A motor failure in 2005 allowed 
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Fig. 1. — Sky maps from the ARCADE 2 2006 flight (thermodynamic temperature in Galactic coordinates). Un-observed regions are 
masked in grey. The Galactic plane is clearly visible. 


only a single channel to view the sky (Singal et al. 2006). 
The 2006 flight successfully obtained observations at 3.3, 
8.3, 10.2, 30, and 90 GHz. Data from the 2006 flight are 
consistent with a blackbody CMB spectrum, but show 
a clear detection of an extragalaetie radio background. 
This paper describes the analysis of Galactic emission 
from the 2006 flight. Companion papers describe the 
ARCADE 2 instrument (Singal et al. 2011), the calibra- 
tion and sky temperature analysis (Fixsen et al. 2011) 
and the results from fits to potential distortions from a 
CMB blackbody spectrum (Seiffert et al. 2011). 

2. OBSERVATIONS AND SKY MAPS 

The ARCADE 2 instrument launched from the 
Columbia. Scientific Balloon Facility in Palestine, TX car- 
rying a complement of 7 Dicke-switched radiometers and 
1800 liters of liquid helium. During ascent, the Dicke 
switch failed in one radiometer. The remaining radiome- 
ters began sky observations on 2006 July 22 at 05:08 UT, 
continuing through 08:11 UT just prior to flight termi- 
nation. 

Figure 1 shows the absolute sky temperatures binned 
by Galactic coordinates. Each radiometer views the sky 
through a corrugated conical feed horn, scaled with wave- 
length to produce a Gaussian beam with 11?6 full width 


at half maximum (Singal et al. 2005). The beams point 
30° from the zenith and the entire payload spins at 0.6 
RPM so that the beams scan a circle of 60° diameter cen- 
tered oil the zenith. The spin frequency is comparable 
to the 1 // knee of the instrument noise, so that strip- 
ing at the mK level is visible in the binned data. Each 
radiometer observed roughly 7 % of the sky. 

The measured sky temperatures are dominated by the 
CMB monopole. Galactic emission is clearly visible in 
the data at 3, 8, and 10 GHz. At higher frequencies, 
the smaller Galactic signal and higher instrument noise 


TABLE i 

Frequency Bands For Foreground Analysis 


Frequency Bandwidth White Noise 
(GHz) (MHz) mK s 1 / 2 


3.20 

210 

11.8 

3.41 

220 

10.1 

7.98 

350 

5.5 

8.33 

350 

5.2 

9.72 

860 

3.0 

10.49 

860 

3.0 



ARCADE 2 Observations of Galactic Radio Emission 


3 


combine to prevent a clear detection of Galactic emission. 
Throughout this paper, we consider only data from the 3, 
8, and 10 GHz radiometers. Table 1 summarizes the in- 
flight performance of each radiometer. We compute the 
instrument noise by comparing the variance of the data 
within each sky bin to the number of observations in that 
bin. This value for the noise thus includes calibration as 
well as the effect of low-frequency striping. 

The observed emission is a combination of CMB, syn- 
chrotron, and free-free emission. Synchrotron emission 
results from the acceleration of cosmic ray electrons in 
the Galactic magnetic field. For a power-law distribu- 
tion of electron energies N(E) oc E~ v propagating in a 
uniform magnetic field, the synchrotron emission is also 
a power law, 

T A (v) ex v~ 0 ‘ (1) 

where T A is antenna temperature, v is the radiation fre- 
quency, and fi 3 = — (p+3)/2 (Rybicki & Lightman 1979). 
Free-free emission is also a power law in antenna temper- 
ature, with spectral index fig = —2.15. 

We evaluate the observed spectral index in the AR- 
CADE 2 sky maps by plotting the temperature observed 
at one frequency against the temperature of the same sky 
pixel observed at a different frequency. Since the AR- 
CADE 2 absolute calibration is derived from comparison 
of sky data to a blackbody calibrator of known physical 
temperature, the sky solutions have units of thermody- 
namic temperature (Fixsen et al. 2011). In order to 
evaluate the spectral index, we convert the sky maps 
(Fig. 1) from thermodynamic to antenna temperature 


exp(.r) — 1 


T 


( 2 ) 


where T is thermodynamic temperature, x = hv/kT, 
h is Planck’s constant, and k is Boltzmann’s constant. 
Figure 2 shows the resulting TT plots. 

The same Galactic features are observed at each fre- 
quency. A linear fit to each plot yields spectral index 
0 = -2.43 ± 0.03 from 3.3 to 8.3 GHz, 0 = -2.47 ± 0.02 
from 3.3 to 10.3 GHz, and fi = —3.12 ± 0.17 from 8.3 
to 10.3 GHz, where the quoted errors include the mea- 
surement uncertainties in both maps for each plot. The 
measured index is consistent with a superposition of syn- 
chrotron and free-free emission. 


3. SPATIAL STRUCTURE 

A search for distortions in the CMB spectrum requires 
correction of Galactic foreground emission. The AR- 
CADE 2 data has only six frequency channels contain- 
ing significant Galactic emission (Table 1), and can not 
uniquely separate the observed emission into individual 
components (synchrotron, free-free, or other sources). 
We combine the ARCADE 2 data with other sky sur- 
veys in order to model foreground emission. 

A widely-used technique describes the observed emis- 
sion at each frequency v as a linear combination of fixed 
“template” maps, 

Ta(v,P ) = y' ti a i (L')Gi(p) (3) 

i 

where p is a pixel index and Gi(p) are a set of sky maps 
tracing different components of the interstellar medium. 
The spectral dependence of the fitted coefficients a(v) 
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Fig. 2. — TT plots from the ARCADE 2 sky maps in units of 
antenna temperature. The best-fit spectral index (solid line) is con- 
sistent with a superposition of synchrotron and free-free emission. 


can then be ascribed to the component traced by each 
template. 

Several choices for template maps tracing free-free 
emission are available. Gold et al. (2009) use a 
maximum-entropy algorithm to derive a map of free-free 
emission based on the Wilkinson Microwave Anisotropy 
Probe (WMAP) 5-year data and a compilation of Ha 
emission (Finkbeiner 2003). Since the maximum- 
entropy method (MEM) template is heavily weighted by 
the microwave data in the Galactic plane, it minimizes 
extinction effects in the Ha map. As alternatives to the 
maximum-entropy template, we also use either a map 
of thermal dust emission (Finkbeiner, Davis, & Schlegel 
(1999) model 8) or the FIRAS map of Gil emission at 
158 //,rn (Fixsen, Bennett, & Mather 1999), neither of 
which is affected by extinction. 

The choice of synchrotron template depends on the 
treatment of spatial variation in the synchrotron spectral 
index. To the extent that this variation can be neglected. 
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TABLE 2 

Coefficients for Spatial Templates 


Frequency 

(GHz) 

Off 

(raK nW -! m 2 sr) 

OL S 

(m K/K) 

3.20 

3.22 ±0.11 

2.02 ± 0.05 

3.41 

3.06 ± 0.09 

1.70 ± 0.04 

7.98 

0.26 ± 0.05 

0.24 ± 0.02 

8.33 

0.28 ± 0.05 

0.24 ± 0.02 

9.72 

0.39 ± 0.03 

0.04 ± 0.01 

10.49 

0.37 ± 0.03 

0.05 ± 0.01 


the 408 MHz survey (Haslain et al, 1981) provides a 
high signal- to- noise ratio tracer for synchrotron emission. 
Since the 408 MHz survey contains both free-free and 
synchrotron emission, we correct it for free-free emission 
using the MEM free-free template scaled to 408 MHz 
using a spectral index /% = —2.15 so that to first order 
the template represents only synchrotron emission. We 
convolve both the synchrotron and free-free templates to 
match the ARCADE 1 1?6 beam width. 

The ARCADE 2 data are dominated by the CMB 
monopole, even towards the Galactic plane. Both the 
Galactic and extragalactic signals also contain monopole 
terms, as do the synchrotron and free-free templates. We 
thus include a monopole template as a nuisance parame- 
ter and compute the template coefficients Q-i(v) by min- 
imizing 


,2 _ ST' \ t a (p-v) - ftp (id - aff(v)Gg(p) - a s (v)(L (p' ~ 

' 2 s a U ,/}2 

p j 

( 4 ) 

where op, an, and a s are the coefficients for the 
monopole, free-free, and synchrotron templates, respec- 
tively. and a(p, v) is the instrument noise in each pixel 
and frequency channel. 

The spatial structure in the ARCADE 2 data may be 
described as a superposition of emission traced by the 
Haslain 408 MHz survey and the FIRAS map of Cn emis- 
sion. Table 2 shows the fitted coefficients for the best-fit 
template combination. Fitting the coefficients from each 
template to a power law in frequency yields spectral in- 
dices —2.5 ±0.1 for emission traced by the 408 MHz sur- 
vey and — 2.0 ± 0.1 for emission traced by the Cl I map. 

The best template model uses the 408 MHz survey to 
trace synchrotron emission and the FIRAS map of Cll 
line emission to trace free-free emission. We have re- 
peated the template fits using different tracers for syn- 
chrotron or free-free emission, obtaining broadly similar 
results. In all cases, the best choice for free- free tem- 
plate is either the Cll map or the thermal dust map. We 
note that spatial structure in the ARCADE 2 maps is 
dominated by t.he Galactic plane, where extinction cor- 
rections for free-free maps derived from Ha emission are 
worst. The ratio of emission traced by the Cll free-free 
template to emission traced by the 408 MHz synchrotron 
template is (Tg/T syn c h) = 0.16 ± 0.01 evaluated for lati- 
tudes |6j < 10° in the lowest ARCADE 2 channel at 3.15 
GHz. The ratio of free-free emission to the total Galactic 
plane emission is (Tff /T lotBL \) — 0.10 ± 0.01. 


4. GALACTIC POLAR CAP TEMPERATURE 

The sky temperature measured by ARCADE 2 in- 
cludes contributions from Galactic, extragalactic, and 
cosmic sources. The 2.7 K cosmic microwave back- 
ground dominates the measured temperatures. Galac- 
tic emission is dominated by synchrotron emission, with 
a smaller contribution from free-free sources. The inte- 
grated contribution of similar synchrotron and free-free 
emission in external galaxies constitutes an isotropic ex- 
tragalactic radio background. Accurate measurement of 
the extragalactic or cosmic backgrounds requires reliable 
determination of both the spatial structure and the zero 
level of the Galactic emission. 

Although the cosmic microwave background may be 
distinguished from Galactic or extragalactic radio emis- 
sion by tiie different frequency dependences, spectral fit- 
ting alone can not distinguish Galactic emission from 
an extragalactic component of similar spectral behav- 
ior. The template model above reproduces the spatial 
structure of the observed Galactic emission, but is in- 
sensitive to emission described by an additive constant 
in each map. If the total Galactic brightness were de- 
termined along some fiducial line of sight, we could add 
a monopole term to the template model to match the 
total Galactic brightness along that line of sight. The 
adjusted model would then characterize Galactic emis- 
sion along all other lines of sight. 



Fig. 3. — Temperature of the 408 MHz survey (top) and AR- 
CADE 3.15 GHz channel (bottom) binned by csc(b) for latitudes 
b > 10° (Northern hemisphere). Sky temperatures include the 
CMB and any extragalactic background. The scatter about the 
best-fit line results from the higher-order spatial structure in the 
maps. For clarity, the plotted uncertainties have been inflated by 
a factor of 20 (408 MHz survey) or 5 (ARCADE); the actual sta- 
tistical uncertainties for the mean of each bin are smaller than the 
symbol size. The data are consistent with a plane-parallel model 
whose measured slope provides an estimate of the Galactic emission 
at the poles. 
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The north and south polar caps (latitude |6| > 75°) 
provide convenient reference lines of sight. A commonly 
used method to estimate total Galactic brightness toward 
the poles models the Galaxy as a simple plane-parallel 
structure and fits the spatial distribution of the radio 
continuum at Galactic latitude |6[ > 10° to the form 

T A (v,p) =c{v)+T g {v) esc \b\ (5) 

where b is the Galactic latitude of each pixel p and 
the constant c(v) accounts for the contribution from the 
CMB or extragalactic backgrounds at each frequency v. 
Figure 3 shows typical results for the northern hemi- 
sphere, both for tiie lowest frequency ARCADE 2 chan- 
nel and the 408 MHz survey. The high-latitude structure 
of the radio sky is consistent with a plane-parallel model, 
whose slope defines the antenna temperature Tq{ v) of 
Galactic emission toward the Galactic poles. 

The simple cosecant, model traces radio emission as- 
sociated with the plane-parallel structure of the Galaxy, 
but averages over significant longitudinal structure and 
is insensitive to emission from a halo or similar extended 
high-latitude structure. As a cross-check, we correlate 
maps of the radio sky with maps of line emission origi- 
nating within the Galaxy to estimate the radio intensity 
corresponding to the measured Galactic line brightness 
along high latitude lines of sight. Several lines could be 
used to trace Galactic microwave emission. Ha emission 
from 3-2 transition in neutral atomic hydrogen has been 
mapped over the full sky and is a well-established tracer 
of free-free emission. However, Ha emission predomi- 
nantly traces free-free emission from the warm ionized 
medium, and may not accurately trace emission from the 
brighter synchrotron component. In addition, Ha emis- 
sion suffers from dust extinction in the Galactic plane 
making it le&s reliable in regions where the Galactic radio 
emission is brightest. The 21 cm fine-structure line from 
neutral hydrogen has also been mapped over the full sky. 
However, since Hi emission originates from the neutral 
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Fig. 4. — Correlation between the ARCADE 3.15 GHz inten- 
sity and the Cn 158 jim atomic line in units nW m -2 sr“*. The 
data clearly divide into two regions each with radio intensity pro- 
portional to the square root of the Cn intensity. The upper track 
originates from pixels at longitude l ~ 45° (closer to the Galactic 
center), while the lower track originates from pixels at longitude 
l ~ 95° (closer to the Cygnus region). 


component of the interstellar medium, it is unlikely to 
trace microwave emission from ionized regions. The ClI 
line at 158 pm wavelength from singly ionized carbon can 
be used to trace diffuse radio emission. It is optically thin 
even in the plane, does not suffer from dust extinction, 
and is an important cooling mechanism for the diffuse 
interstellar medium. It has been mapped over the full 
sky bv the COBE/FIRAS instrument (Fixsen, Bennett, 
& Mather 1999). 

Figure I shows the correlation between the ARCADE 
3.15 GHz map and the Cn 158 pm map. The ClI map 
has been smoothed to angular resolution 11?6 to match 
the ARCADE resolution. The data segregate into two 
tracks, an upper track from emission near the Galac- 
tic center and a lower track from emission closer to the 
Cygnus region. Both tracks show radio emission propor- 
tional to the square root of the ClI intensity. This would 
be expected if the radio emission were dominated by syn- 
chrotron emission (proportional to the density n in the 
interstellar medium), but with collisionally excited ClI 
intensity proportional to n 2 . Provided both the high- 
latitude radio continuum and the ClI line intensity pre- 
dominantly originate from a small number of emitting 
regions along each line of sight, the radio temperature 
T and line intensity I would then scale asToc/* with 
K = 1/2. Emission along the same line of sight from 
a large number of regions with different local densities 
would force the relation closer to linear (k ■— 1), as would 
a significant admixture of radio free-free emission (with 
antenna temperature proportional to n 2 ) The observed 
dependence shows n in the range 0.5 - 0.7, as can be 
verified by the lack of significant curvature in Figure 4, 
The absence of any such curvature provides additional 
evidence that, free-free emission is faint compared to syn- 
chrotron at the ARCADE 2 bands. Note, though, that 
if a separate template is provided for the synchrotron 
component, the Cl I map may then be used to trace the 
fainter free-free component (at least at the ARCADE an- 
gular resolution). 

The prominent segregation in Fig 4 results at least in 
part from the limited sky coverage. The brightest re- 
gions in the ARCADE 2 maps occur at the portion of 
the ARCADE 2 sky coverage closest to the Galactic cen- 
ter region and again at the portion closest to the Cygnus 
region. The ratio of Cn intensity to the far-infrared con- 
tinuum is known to differ by a factor of two between 
these regions (Fixsen, Bennett, & Mather 1999). The 
radio/Cli correlation for the ARCADE 2 sky coverage 
is dominated by two segregated tracks. Fitting a single 
slope to the ARCADE 2 data could induce a systematic 
error as small shifts in sky coverage view more of one 
track than the other. We minimize this effect by fitting 
two components to the radio/Cli correlation plots, 

2 

T A (v,p) = Y,bi{v) + a,{u){ Icip) ) 0 ' 5 , (6) 

i—1 

where T A is the antenna temperature at frequency v. Ic 
is the Cn intensity, and p denotes a pixel index. After a 
first guess at the fit parameters (intercept and slope for 
each line), we assign each point to one line or the other 
based on the difference in antenna temperature between 
that point and either line. We then re-compute the fit 
parameters for each line using only those points assigned 
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to that line, and then re-assign all points based on the 
new fit parameters. Several iterations suffice to produce a 
stable solution independent of the initial guess. We then 
define the radio/ Gil correlation slope as the arithmetic 
mean of the two fitted components, (a(i/)) = 0.5[oi(i/) + 
a- 2 (V)], with uncertainty broad enough to bracket the two 
cases. The mean value thus derived is consistent with 
similar analysis fitting a single slope to the full sky but 
is less sensitive to small shifts in the sky coverage. Table 
3 shows the radio/Cri correlation. 

The observed correlation between radio emission and 
ClI intensity allows an independent determination of the 
associated Galactic radio emission toward the Galactic 
poles, 

Tg(u) = (oM) (/ cap f, (7) 

where J cap is the mean intensity of the Cii map at the po- 
lar caps (Galactic latitude |6| > 75°) and we take « = 0.5. 
Higher values of k within the observed range 0.5 - 0.7 
lower the estimated Galactic emission by less than the 
half the model uncertainty. Since the correlation coef- 
ficient a(i/) is derived using spatial fluctuations in the 
observed radio and line intensity, this method accounts 
at least roughly for both longitudinal as well as latitudi- 
nal structure on the sky. However, to the extent that the 
radio emission is not co-located with the line emission, 
the correlated emission does not uniquely determine the 
zero level of either component. 

We repeat both the cosecant fit (Eq. 5) and the 
radio/Cn correlation (Eq. 7) independently for the 
northern and southern hemispheres, using each of the 
ARCADE 2 sky maps as well as a selection of lower- 
frequency radio surveys. Surveys at 22 MHz (Roger et 
al. 1999), 45 MHz (Maeda et al. 1999; Alvarez et al. 
1997), 408 MHz (Haslam et al. 1981), and 1420 MHz 
(Reich, Testori, & Reich 2001; Reich k Reich 1986) 
have full or nearly full sky coverage at frequencies where 
the sky brightness is dominated by Galactic radio emis- 
sion. For each sky map, we compute the cosecant slope 
Tg(^) and its uncertainty, including both statistical and 
calibration uncertainties. Although the binned data are 
broadly consistent with a plane-parallel model, higher- 
order structure causes the binned data to show greater 
scatter about the best-fit slope than would be expected 
given the formal statistical uncertainty in the mean of 
each bin. We account for higher-order structure by in- 
flating the statistical uncertainty to force \' 2 to unity 
per degree of freedom about, the best-fit. line. We then 
add the calibration uncertainty (of order 10% for the 
low-frequency radio surveys) in quadrature with the (in- 
flated) statistical uncertainty to derive the total uncer- 
tainty in the cosecant slope for each sky survevy. 

Figure 5 compares the Galactic temperature towards 
the north polar cap, derived from mean Cn correlation, 
to the value derived from the cosecant fit. Roth methods 
provide similar estimates for the Galactic component of 
the total brightness temperature towards the north polar 
cap. This Galactic component is consistent with a single 
power law 

Tg(is) = Tg&,(i j /v o } 0 (8) 

with spectral index 8 — —2.55 ± 0.03 and ampli- 
tude To,,! = 10.12 ± 0.90 K at reference frequency 



Fig. 5. — Galactic emission Tq(v) towards the North galactic 
pole, derived from a cosecant fit to the spatial structure (open 
circles) and the radio/Cn correlation (filled circles). Both methods 
agree and are consistent with a single power law over the frequency 
range 22 MHz to 10 GHz. The fitted Galactic temperature towards 
the polar caps can be combined with the template model of the 
spatial structure to fully specify the Galactic model. 


t'o = 0.31 GHz 9 . 

Several mechanisms could alter the total high- latitude 
Galactic emission from a power law form. Energy losses 
from cosmic rays propagating in the interstellar medium 
will cause a progressive steepening of the synchrotron 
spectrum at higher frequencies. Comparison of the 408 
MHz survey with data from the Wilkinson Microwave 
Anisotropy Probe (WMAP) mission show some evidence 
for such steepening at frequencies above 20 GHz (Gold et 
al. 2009; Bennett et al. 2003). Free-free absorption can 
become appreciable at sufficiently low frequencies, lower- 
ing the total temperature to cause a shallower spectrum 
below ~50 MHz where the synchrotron antenna temper- 
ature exceeds the typical 8000 K electron temperature in 
the diffuse interstellar medium. We may use the modeled 
Galactic temperatures to examine such effects. 

The data in Figure 5 are not independent. Data from 
different ARCADE bands share common-mode ther- 
mometry uncertainty and small corrections for instru- 
mental emission. The gain calibration of the 22 MHz 
survey is derived in part, from comparison of the 22 MHz 
survey to the 408 MHz survey, thereby correlating these 
results. We account for these correlated uncertainties 
when deriving the fitted parameters. Including these cor- 
relations, the distribution of the individual points in Fig. 
5 with respect to the best-fit power- law model is con- 
sistent with the measurement uncertainties. Additional 
free parameters (e.g. for possible spectral steepening) 
are not statistically justified. 

The analysis above uses a single map (Co line emis- 
sion) to trace Galactic emission from the diffuse inter- 
stellar medium. If a significant fraction of Galactic radio 
emission originated from a component of the interstellar 
medium not well sampled by ClI emission, the result- 
ing correlation w'ouid under-estimate the total Galactic 

9 We chose reference frequency 0.31 GHz to minimize the co- 
variance between the fitted amplitude and spectral index. See the 
discussion in Pixsen et a!. (2011). 
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TABLE 3 

Radio/Cii Court elation Slope 

Frequency Correlation Slope a (inK [ nW rn 2 sr 'M'" 0,5 j 


(GHz) 

Component, 1 

Component, 2 

Mean 


0.022 

(1.12 ±0.06) x 

: 10 4 

(0.74 ± 0.02) x 

: 10 4 

(0.93 ±0.19) x 

ilcF' 

0.045 

(2.74 ± 0.07) x 

: 10 3 

(1.56 ±0.03) x 

: 10 3 

(2.15 ±0.63) X 

: 10 3 

0.408 

18.1 ±0.3 



8.5 ±0.1 


13.3 ±5.0 


1.420 

(5.11 ±0.08) x 

10-1 

(2.94 

± 0.04) X 

io- 1 

(4.02 ± 1.10) x 

IO” 1 

3.20 

(5.63 ±0.12) x 

10 2 

(2.88 

± 0.09) x 

10~ 2 

(4.25 ± 1.38) x 

10“ 2 

3.41 

(5.05 ±0.11) x 

10-2 

(2.61 

± 0.08) x 

io- 2 

(3.82 ± 1.20) x 

10- 2 

7.98 

(5.73 dt 0.37) x 

10- 3 

(1.99 

± 0.36) x 

10- 3 

(3.86 ± 1.87) x 

10-3 

8.33 

(6.32 ± 0.34) x 

lO- 3 

(2.59 

zb 0.26) x 

10“ 3 

(4.46 ± 1.86) x 

10-3 

9.72 

(3.68 ±0.15} x 

io- 3 

(2.37 

±0.16) x 

10-3 

(3.03 ± 0.66) X 

10 “3 

10.49 

(3.67 ±0.16) x 

10 3 

(2.17 

±0.18) x 

10-3 

(2.92 ±0.75) X 

10~ 3 


emission toward the polar cap. We test for additional mi- 
crowave emission from other components of the interstel- 
lar medium by repeating the correlation analysis using a 
simultaneous fit to three line maps chosen to sample dif- 
ferent components of the diffuse interstellar medium: the 
(square root of the) COBE/FIRAS map of Cll emission, 
a map of Ha emission (Finkbeiner 2003), and the Lei- 
den/ Argentine /Bonn Galactic Hi Survey (Kalherla et al. 
2005). The simultaneous fit produces results nearly iden- 
tical to a fit using just the Cll map. shifting the derived 
Galactic polar cap temperature by 0.2 ± 0.4 K at refer- 
ence frequency 0.31 GHz. Since the limited ARCADE 2 
sky coverage increases the effect of covariance between 
the different line maps, we use the results derived from 
the single Cl I correlation. 

4.1. Composite Galactic Model 

Removing Galactic emission from the ARCADE 2 data 
requires a model for both the spatial structure and the 
total Galactic emission. At each ARCADE frequency, wc 
generate a full-sky model of the spatial structure using 
full-sky template maps multipled by the coefficients fit: ted 
within the region observed by ARCADE 2 (Table 2). We 
then add a constant to the template model to match the 
total Galactic emission derived from the Cll and esc b\ 
fits towards selected reference positions. 

We independently model the total Galactic emission 
toward the north and south Galactic polar caps (|6| > 
75°). For each reference position, we compute the to- 
tal Galactic emission associated with either the plane- 
parallel structure (Eq. 5) or the Cll emission (Eq. 7). 
We apply each method to the ARCADE 2 and radio sur- 
veys, and parameterize the resulting multi-frequency re- 
sults using a power-law model (Eq. 8). As a cross-check, 
we also estimate the total Galactic emission for the cold- 
est, patch in the Northern hemisphere, consisting of all 
pixels within 15° of b = 48°, l — 196° (north of the Galac- 
tic anti-center). Since by definition a mid- latitude cold 
spot is inconsistent with a plane-parallel structure, the 
estimate for this position is based only on the radio/ClI 
correlation. Table 4 shows the estimated Galactic nor- 
malization Tg al and spectral index /3 for each reference 
position. We obtain similar results using either the spa- 
tial morphology or the radio/ClI correlation. 

We fix the offset of the template model (Eq. 3) by com- 
puting the temperature of the template model towards 
these same regions, then adding a constant ao(v) to the 


template model at each frequency v to force the template 
model to match the power-law model along these lines of 
sight. The resulting composite model of Galactic emis- 
sion can be used to correct the calibrated time-ordered 
data in order to estimate the CMB monopole tempera- 
ture and the extragalactic radio background (Fixsen et 
al. 2011). The offsets derived from the three indepen- 
dent regions agree within 5 mK at 3 GHz, which we adopt 
as the uncertainty in the zero level of composite Galactic 
model. 

The spectral index derived for the Galactic emission to- 
ward the polar caps or coldest region, j3 as —2.57 ± 0.03, 
is consistent with the spectral index (3 — —2.5 ± 0.1 de- 
rived from the synchrotron template fit only to ARCADE 
2 data, indicating that high-latitude Galactic emission is 
dominated by synchrotron emission at frequencies below 
10 GHz. A synchrotron spectral index of -2.57 is signifi- 
cantly flatter than the values derived from measurements 
at higher frequencies, but is consistent with the majority 
of measurements below 10 GHz (see, e.g., the tabulation 
in Rogers k. Bowman (2008)). Notable exceptions are 
the measurements by Tartari et, al. (2008) and Platania 
et al. (1998), which prefer 3 « —2.8 for the synchrotron 
component. Both of these results are derived from a 
single strip at declination 5 = +42° and may not be 
representative of larger regions of the sky. 

Although a plane-parallel model provides a good de- 
scription of the latitudinal dependence of the radio sky, 
the projected polar cap temperatures from a esc |6| model 
are well below the actual measured values. For example, 
the North Galactic pole has measured temperature 19±3 
K at 408 MHz, while a esc |6| fit. to the same 408 MHz 
map predicts a polar contribution of only 5.1±0.6 K. Sim- 
ilar results apply to the Southern hemisphere, where the 
measured polar cap temperature of 21 ±3 K significantly 
exceeds the value 4.0 ± 0.5 K obtained from a esc \b\ fit. 
Only 2.7 K of the difference can be attributed to emission 
from the cosmic microwave background, leaving a large 
(12 K) residual. Several authors have noted the existence 
of this bright, isotropic component of high-latitude radio 
emission (Phillips et al. 1981; Beuermann, Kanblach, 
& Berkhuijsen 1985; Sun et al. 2008; Everett, Schiller, 
& Zweibel 2010). The Galactic model presented here 
is consistent, with these models. The high latitude sky 
is dominated be an emission component not associated 
with the plane-parallel structure of the Galaxy, but resid- 
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TABLE 4 

Galactic Emission® Along Selected Lines of Sight 


Parameter Technique 

North Polar Gap 

South Polar Cap 

Coldest Patch 

Radio/Cn 

9.69 ±2.14 

8.83 ±2.07 

5.40 ± 3.86 

T Gal (K) esc |6[ 

10.21 ± 0.99 

8.29 ±0.90 

~~ 

Weighted Mean 

10.12 ±0.90 

8.38 ±0.83 

5.40 ± 3.86 

Radio/Cll 

-2.52 ± 0.07 

-2.53 ± 0.07 

-2.52 ±0.11 

0 esc |6i 

-2.55 ±0.04 

-2.64 ± 0.05 

— 

Weighted Mean 

-2.55 ± 0.03 

-2.63 ± 0.03 

-2.52 ±0.11 


Galactic emission Tq(u) = 7’ofti(i//Vo}^ with reference frequency uq = 0.31 GHz, 


ing in either a Galactic halo or an isotropic extragalactic 
background. 

5. SPINNING DUST 

The ARCADE 2 data can also be used to search for ad- 
ditional microwave emission associated with interstellar 
dust. The microwave sky is known to contain a compo- 
nent spatially correlated with far-infrared dust emission 
but not with synchrotron-dominated surveys at 408 or 
1 120 MHz (Kogut et al. 1996a,b; de Oliveira-Cost, a et 
al. .1997; Leitch et al. 1997; Bennett et al. 2003). The 
frequency spectrum of the correlated component differs 
markedly from the emission spectrum of thermal dust., 
showing a spectral index f3 « —2.2 in antenna tempera- 
ture from 20 to 50 GHz. 

Two main candidates have emerged to explain this 
“anomalous” correlated emission component. Electric 
dipole emission from a population of small, rapidly ro- 
tating dust grains will produce a spectrum with a broad 
peak in the frequency range 20-40 GHz (Draine &: Lazar- 
ian 1998). Miville-Deschenes et al. (2008) analyze data 
from the WMAP and conclude that spinning dust ac- 
counts for the majority of emission in the Galactic plane 
at 23 GHz. An alternative model explains the correlated 
emission as flat-spectrum synchrotron emission associ- 
ated with star formation activity (Bennett et al. 2003; 
Gold et al. 2009). Energy losses as electrons propagate 
from their origin steepen the synchrotron index away 
from these sources. Since a flat-spectrum component will 
increasingly dominate at higher frequencies, the spatial 
morphology of synchrotron emission at higher frequen- 
cies should increasingly resemble thermal dust emission 
from the same star formation activity. 

Efforts to distinguish between these models have 
largely focused on the frequency spectrum of the dust- 
correlated component inferred by correlating a map trac- 
ing thermal dust emission against microwave maps at dif- 
ferent observing frequencies. The correlation with ther- 
mal dust emission is seen to peak at frequencies between 
20 and 30 GHz, falling at frequencies below 23 GHz 
(de Oliveira-Costa et al. 2004; Fernandez- Cerezo et al. 
2006; Hildebrandt et al. 2007). This falling spectrum 
at low frequencies has been cited as evidence favoring 
spinning dust models. 

The decrease in dust-correlated emission at frequencies 
below 20 GHz, although consistent with spinning dust 
emission, does not necessarily support such models over 
the flat-spectrum synchrotron alternative. Both models 
predict weaker correlation with thermal dust emission at 
frequencies below 20 GHz. For the spinning dust model. 


the weaker correlation results from the lower amplitude 
of the spinning dust emission. The flat-spectrum syn- 
chrotron model, however, also predicts a weaker correla- 
tion at lower frequencies, since by construction the spa- 
tial morphology of the sky in this model varies smoothly 
from the WMAP data at 23 GHz (observed to correlate 
well with the thermal dust- morphology) to the Haslam 
map at 408 MHz (observed to correlate only weakly with 
thermal dust emission). Model discrimination based on 
cross-correlations with thermal dust emission must rely 
on rigorous statistical comparison of competing models, 
and not simply on the general trend toward weaker cor- 
relations at lower frequencies. 

A simpler test relies on the absolute spectrum of Galac- 
tic emission. At frequencies below 30 GHz, thermal dust 
emission is negligible, so that the total Galactic emis- 
sion becomes a superposition of free-free, synchrotron, 
and spinning dust emission. Both free-free and syn- 
chrotron emission increase monotonically at lower fre- 
quencies. Spinning dust, in contrast, decreases in ampli- 
tude below 20 GHz, so that the spectrum of the com- 



Fig. 6. — Synchrotron and spinning dust emission for two models 
with different spinning dust amplitude. The red (blue) curves show 
emission from a model with high (low) dust normalization. Dashed 
lines show the dust emission, dotted lines show synchrotron emis- 
sion. and solid lines show the combined emission from each model. 
The amplitude of the combined emission is fixed at both 23 GHz 
and 408 MHz. Measurements below 10 GHz are dominated by 
synchrotron, which can be used to infer the dust normalization; 
models with more dust emission at 23 GHz have lower synchrotron 
emission. Grey bars indicate the ARCADE 2 frequency bands. 
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billed Galactic emission below 20 GHz can be used to 
place limits on the contribution of spinning dust without 
reference to detailed spatial correlations with a thermal 
dust, template. 

Figure 6 illustrates the concept. Consider Galactic 
emission consisting of a superposition of synchrotron and 
spinning dust (ignoring for the moment the smaller con- 
tribution from free-free emission). Following Miville- 
Deschenes et al. (2008), we specify the spinning dust 
amplitude as a fraction of the total Galactic brightness at 
23 GHz, assumed here to represent the peak in the dust 
spectrum. As the dust normalization is increased, the 
synchrotron amplitude at 23 GHz must decrease to keep 
the total emission constant. The spinning dust spectrum 
falls rapidly, so that emission below 10 GHz is dominated 
by the synchrotron component. A model with high dust 
normalization (red curves in Fig. 6) will thus have fainter 
synchrotron emission at frequencies of a few GHz, while 
a model with less spinning dust (blue curves) will have 
brighter synchrotron emission. This shift in the ampli- 
tude of Galactic emission at frequencies below 10 GHz 
is a sensitive test for spinning dust, without resort to 
detailed spatial correlations. 

We implement this test using a simple model of Galac- 
tic emission. We assume that the spatial distribution 
of spinning dust emission is traced by a template map 
of thermal dust emission, and fix the amplitude of as- 
sociated spinning dust emission by scaling the template 
map so ( hat the re-scaled dust map forms a specified 
fraction of the Galactic plane intensity measured by the 
WMAP 23 GHz map. This allows a simple normalization 
of spinning dust in terms of its relative contribution to 
the total Galactic plane intensity at 23 GHz. We use the 
COBE/DIRBE 240 pm map of thermal dust emission 
(Reach et al. 1996) as the thermal dust template since 
it is dominated by thermal dust emission but unaffected 
by extinction in the plane. For specificity, we define the 
Galactic plane mask using all pixels lying within the AR - 
CADE 2 observation pattern with latitude \b\ < 20°, and 
use this mask for all computations so that the model re- 
sults are not affected by differing sky coverage between 
ARCADE 2 and other surveys. 

We then model the Galactic emission spectrum as fol- 
lows. We first correct the WMAP 23 GHz map and 
the Haslam 408 MHz map by subtracting the WMAP 
maximum-entropy model of free-free emission (Gold et 
al. 2009) using a spectral index ;% = —2.15. The cor- 
rected maps then contain only synchrotron and (possi- 
bly) spinning dust emission. We assume that the spin- 
ning dust spectrum peaks at 23 GHz and scale the nor- 
malized spinning dust map to lower frequencies using the 
Draine & Lazarian (1998) model for the warm neutral 
medium. After correction for free- free and spinning dust 
emission, the 23 GHz and 408 MHz maps contain only 
synchrotron emission, which we use to define the syn- 
chrotron amplitude and spectral index for each pixel. 

The resulting model (synchrotron, free-free, and spin- 
ning dust) can be used to estimate the combined Galac- 
tic emission at frequencies between 408 MHz and 23 
GHz. We compare the model to the ARCADE 2 data 
by smoothing the model map of combined emission to 
the ARCADE 2 angular resolution, and then computing 
the mean Galactic emission of the smoothed model for 
all pixels within the ARCADE 2 galactic plane mask. 


Figure 7 show's the result. If spinning dust is negligible, 
the synchrotron contribution at 23 GHz is maximal and 
the model approximates the fiat-spectrum synchrotron 
model. As the spinning dust amplitude increases, the 
synchrotron contribution at 23 GHz decreases and the 
synchrotron spectral index steepens. The combination 
of lower synchrotron amplitude and falling dust spec- 
trum combine to lower the total model emission across 
the ARCADE 2 frequency bands. The ARCADE 2 data 
lie below the model prediction for no spinning dust, and 
are consistent with spinning dust contributing 0.4 ±0.1 
of the total K-band Galactic plane emission. 

WMAP is a differential instrument and is insensitive 
to any monopole emission component. The zero level of 
the 23 GHz map is set using the cosecant dependence 
on Galactic latitude (Hinshaw et al. 2009). The 408 
MHz survey and the ARCADE 2 data, however, both 
include monopole contributions. To prevent discrepant 
treatment of the map zero levels from affecting the model 
predictions, we remove a monopole from the 408 MHz 
survey and the ARCADE 2 sky maps, and re-set the 
zero level of each map using a esc |6| fit to the spa- 




Fig. 7. — Mean intensity of Galactic plane emission for the 
ARCADE 2 data, compared to model predictions with and without 
spinning dust (see text). Top panel: Model predictions and data. 
The upper solid curve shows the model prediction for no spinning 
dust, averaged over the ARCADE 2 sky coverage with |fe| < 20 °. 
The lower solid curve shows the model prediction for spinning dust 
amplitude equal to 60% of the total Galactic plane emission at 23 
GHz. The dotted curve show's the best fit to the ARCADE 2 data. 
The spectrum of emission from spinning dust is assumed to peak 
at 23 GHz, so a higher spinning dust fraction at 23 GHz produces 
lower total emission at the ARCADE 2 frequencies. The bottom 
panel shows the same data, normalized by dividing each point by 
the model prediction for no spinning dust. Galactic plane emission 
observed by ARCADE 2 is consistently fainter than expected for 
a model with no spinning dust, and is consistent with spinning 
dust contributing 0.4 ±0.1 of the total Galactic plane emission at 
reference frequency 23 GHz. 
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tial structure in each map to match the processing of 
the WMAP data. We restrict analysis to pixels at low 
Galactic latitude (|6| < 20°) where Galactic emission is 
brightest so that uncertainties in the zero level have min- 
imal effect. The largest uncertainty in the analysis is 
the correction for free-free emission. We use the WMAP 
maximum-entropy model of free-free emission, derived 
assuming that spinning dust contributes negligibly to the 
total Galactic emission. We test the sensitivity of the re- 
sult to the free-free normalization by repeating the anal- 
ysis using the same maximum-entropy model scaled by 
a constant normalization factor. Changing the free-free 
normalization by a factor of 2 affects the best- fit spin- 
ning dust normalization by approximately 0.08, and is 
included in the total quoted uncertainty. 

6. CONCLUSIONS 

We use the ARCADE 2 absolutely calibrated obser- 
vations of the sky to model Galactic emission at fre- 
quencies 3, 8, and 10 GHz and angular resolution 11?6. 
TT plots of the binned sky maps show net spectral in- 
dex /3 = —2.43 ± 0.03 between 3.3 and 8.3 GHz and 
[i — —2.47 ±0.02 between 3.3 and 10.3 GHz, consistent 
with a superposition of synchrotron and free-free emis- 
sion, The spatial structure in the maps can be described 
using two spatial templates, the Haslam 408 MHz sur- 
vey to trace synchrotron emission and the COBE/FIRAS 
map of Cll emission to trace free-free emission. Fitting 
these templates to the ARCADE 2 maps yields spec- 
tral index ,i,ynch — —2.5 i 0.1 for emission traced by 
the synchrotron template, with free-free normalization 
(Tff/Xtotai) — 0.10+0,01 evaluated for latitudes \b\ < 10° 
in the lowest ARCADE 2 channel at 3.15 GHz. 

The template model only specifies Galactic emission 
up to an additive constant. We fully specify the Galac- 
tic model by computing the temperature of the tem- 
plate model toward the north and south polar caps, 
then adding a constant to the template model to match 
the Galactic temperature derived from a simple plane- 
parallel model. We extend the analysis to full-sky sur- 
veys at; lower frequencies and find that Galactic polar cap 
temperature is consistent with a single power law over the 
frequency range 22 MHz to 10 GHz, with spectral index 
f3 — —2.55 ± 0.03 and normalization 10.12 ± 0.90 K at. 
reference frequency vq — 0.31 GHz. Emission associated 
with the plane-parallel structure accounts for only 30% 
of the observed high-latitude sky temperature, with the 


residual in either a Galactic halo or an isotropic extra- 
galactic background. 

ARCADE2 produces well-calibrated maps of Galactic 
emission. We use these maps to test for contributions 
from spinning dust near the Galactic plane. Previous 
tests for spinning dust at frequencies below 20 GHz have 
used only the spatial correlation with template maps 
tracing thermal dust emission. The frequency depen- 
dence of the dust-correlated component is not a strin- 
gent test for spinning dust, since a weaker correlation 
at lower frequencies is expected both for models with 
significant spinning dust contributions as well as mod- 
els with flat-spectrum synchrotron but no spinning dust. 
We use a simple model of the total Galactic emission to 
predict the mean intensity of the Galactic plane in the re- 
gion observed by ARCADE 2. The model normalizes the 
spinning dust contribution at 23 GHz, and computes the 
expected total contribution from free-free, synchrotron, 
and spinning dust as a function of frequency and spinning 
dust normalization. The ARCADE 2 data consistently 
show less emission in the Galactic plane than a model 
with no spinning dust, and are consistent with spinning 
dust contributing 0.4±Q,1 of the Galactic plane emission 
at 23 GHz. 
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